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Abstract 

We investigate the number of symmetric matrices of non-negative integers 
with zero diagonal such that each row sum is the same. Equivalently, these are 
zero diagonal symmetric contingency tables with uniform margins, or loop-free 
regular multigraphs. We determine the asymptotic value of this number as the 
size of the matrix tends to infinity, provided the row sum is large enough. We 
conjecture that our answer is valid for all row sums. 



1 Introduction 



Let M{n, i) be the number of n x n symmetric matrices over {0, 1, 2, ... } with zeros 
on the main diagonal and each row summing to i. 

Define A = i/{n—l), which is the average matrix entry ignoring the diagonal. Our 
interest is in the asymptotic value of M{n, £) as n — )■ oo with i being a function of n. 

Very little seems to be known about this problem. The asymptotic value of M{n, 3) 
was determined by Read in 1958 [13]. According to Bender and Canfield [3], de Bruijn 
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Figure 1: An example of a matrix counted by M(9, 20) = 1955487489759152410696. 



extended this to M(n, i) for fixed £ but failed to publish it. In any case, |3] generalised 
the result to bounded but possibly non-equal row sums. The method of [7] applied to 
the intermediate results of [12] is likely to reveal the asymptotic value of M{n, t) when 
£ = o(n^/^), but this calculation is yet to be performed. 

In this paper we treat the case of large C. and manage to find the asymptotics 
whenever £ > Cn/\ogn for any C > |. We will use the multidimensional saddle 
method, which was previously applied successfully to the corresponding {0, 1} problem 
by McKay and Wormald pjj and to the corresponding non-symmetric problem by 
Canfield and McKay [5]. For the non-symmetric problem with mixed row and column 
sums, see Barvinok and Hartigan [1]. 

Our theorem is as follows. 

Theorem 1. Let a and b be positive real numbers such that a + b < |. Let i = i{n) be 
such that in is even and X = i/{n — 1) satisfies 



3a log n 

Then as n ^ oo, 



In Section [21 we express M{n,i) as an integral in n-dimensional complex space and 
divide the domain of integration into three parts, then in Section [3] we estimate the 
integral in two of the parts. In Section HJ we show that the third part is negligible 
in comparison, for the case that i is bounded by a polynomial in n. In Section O we 
complete the proof for large i by means of the theory of Ehrhart quasipolynomials. 

In Section [6], we show that the value given in Theorem [1] closely matches that 
given by a naive probabilistic model, and conjecture with numerical evidence that this 
relationship holds for all i > 0. 
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Throughout the paper, asymptotic notation hke 0{f{n)) refers to the passage of n 
to oo. We will also use a modified notation 0{f{n)). A function g{n) belongs to this 
class provided that 

g{n)=Oif{n)n'^'), 
where the c is a numerical constant (independent of n and e). 



2 An integral for M(n, i) 



We now express M{n, t) as an integral in ra- dimensional complex space and outline a 
plan for estimating it. 

We begin with a generating function in n variables xi, . . . , Xn, 

n {'^-xjXk)~\ 

l<j<k<n 

for which the coefficient of x{^ ■ ■ ■ x^" is the number of n x n symmetric matrices over 
{0, 1,2,...} with zeros on the main diagonal and row sums £i, . . . , £„. In particular, 
M(n, i) is the coefficient of ■ ■ ■ 

Applying Cauchy's Integral Formula we have 



where each variable is integrated along a contour circling the origin once in the anti- 
clockwise direction. It will suffice to take the contours to be circles; specifically, we will 
put Xj = re*^^ for each j, where 



1 + A 

The reason for our choice of r will become clear in Section [3l 
This gives 

where 

I{n)= ■■■/ - aj:j_,e, (3) 

Let F{0) be the integrand in ([3]). 

The quantity (1 — A(e*(^J+^'=) — 1))^^ has greatest magnitude when 9j + 9k E {0, 2%}, 
and hence F{6) has greatest magnitude when either 6j + 6^ = or 6j + 0^ = 2tt for 



3 



each pair j, k. It is trivial to show that these equations have only two solutions: 6j = 
for all j, and 6j = vr for all j. We will show that the value of J(n) comes mostly 
from the neighbourhoods of these two points; specifically, it comes from two boxes 
T^o,T^n ^ [—^jTr]" defined, for sufficiently small 

no = {6: \ej\ < + X)-^ for all j }, and 

= {6» : lOj + 7r| < + X)-^ for all j }, 

where 6j + vr is taken mod 2tt. Note that the operation 6j 6j + n for all j, which 
maps TZq and TZj^ onto each other, preserves F{6) since ni is even. Also note that 
TZq n TZj, = 0. We denote the region outside of the boxes as 

7^^ = [-7^,7^]"\(7^oU7^^). (4) 

If X C [-vr, vr]", then we let Ix{n) = F{9) dO. For A = O(n^) we will evaluate 
the integral I{n) defined in (|3]) in the following way: 

I{n) = Ino{n) + InA'") + In-{n) 
= 2In,{n) + 0{l) [ \F{e)\dO 

= 2In'{n) + 0{l) f \F{e)\de (5) 
for any 7^' with 7^o C 7^' C [-tt, tt]" \ 7^^. 



3 The main part of the integral 

In this section we estimate the value of the integral J(n) in a convenient region TZ' that 
contains TZq. We begin by quoting several results required for the calculation. 

The following theorem, simplified from |9j, estimates the value of a certain multi- 
dimensional integral. 

Theorem 2. Let e' , e" ,e"' ,e,e, A be constants such that < e' < e" < e'" , e > 0, 
e > 0, and < A < 1. The following is true if e'" and e are sufficiently small. 

Let A = A{n) be a real-valued function such that A{n) = Q{n^'^'). Let B = B{n), 
C = C{n), E = E{n), F = F{n), G = G{n), H = H{n), and I = i{n) be complex- 
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valued functions of n such that B, C, E, F, G,H,i = 0{rf) . Suppose that 

(n n n n 

—An ^ z^j + Bn ^ z^^ + C ^ Zjzl + DrT^ ^ ZjZkZp 
j=l j=l j,k=l j,k,p=l 

n n n n 

j=i j)k=i i)fc=i jM,p=i 



+ in~^/^ ^ ZjZkZpZq + 6{z)\ 



j,k,p,q=l 

is integrable for z = {zi, Z2, ■ ■ ■ , Zn) € [/„ and 5{n) = maXzeUn\^i^)\ = where 
Un = {zC R" : I^jI < n"^/^+^" /or 1 < j < n}, 

where i = e{n) satisfies e" < 2e < e'" . Then, provided the 0{) term in the following 
converges to zero, 

J f{z) dz = l^-^j exp (q^ + O {n-'/^+' + {n^^ + 6{n))Z) 



where 



15B^ SBC 3E F 

Qi = — — + + — - + ^ + 



16^3 8^3 16^3 4^2 4^2 

^_ fl5lm{Bf + Qlm{B)lm{C) + lm{Cf 
V 16^3 



The following lemma defines a linear transformation, adapted from 
Lemma 3. Define c and z = {zi, Z2, . . . , Zn) by 



n 

C 



{l + X)e, = z,--J2^k {l<J<n). (7) 
fe=i 

The transformation 6 = T{z) defined by ([7]) has determinant (1 — c)/(l + A)". For 
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m 



> 1, define fim = Yl^=i^T' ^^^'^ have the following translations. 

n 

(l + A)5^e, = (l-c)/xi 



"■1' 



(1 + \f (^J- + ^^)' = - 2)/i2, 

l<j<k<n 

(1 + A)3 J2 i^j + = in- 4)/i3 + (3(1 - 2c) + 12c/n)fi,fi2 

^-^^'"-'^ + ((_6c + I2c^ - Ac^)/n - 40^(3 - c)/n^)fj,l, 

(1 + A)^ J2 (^J- + ^fe)^ = (^~ 8)/i4 + 3/i^ + (4(1 - 2c) + 32c/n)/ii/i3 

- (24c(l - c)/n + 4:8c^/n'^)f4fX2 
+ (8c2(l - c)(3 - c)/n'^ + 8c3(4 - c)/n^)fii 

From Taylor's Theorem with remainder we have 
Lemma 4. For all real X , 

(1 - A(e^^ - 1))"^ =exp(AiX - |A(1 + \)X'^ - \i\(l + A)(l + 2A)X=' 

+ ^A(l + A)(l + 6A + 6A2)X^ + 0((A + A^)X5)). 

We now present the main result of this section. 

Theorem 5. Under the conditions of Theorem [3, there is a region TV such that 
7^0 ^ 7^' C 37^o C [-71, Tif \ and 

T ( ^ 1 f V^' / 14A^ + 14A-1 , 



v/2VA(l + A)n; 12A(A + 1) 

Proof. Consider the transformation 6 = T{z) defined by ([7]). Define 
n^ = {z: \zj\ < 2n-i/2+"} and 7^' = ^(7^^). 

From ([7]) we have 

\9j\ < y for all j =^ \zj\ < (1 + A)(l - c)-^y for all j, 
\zj\ < y for all j =^ \9j\ < (1 + A)"^(l + c)y for all j. 

These imply, for n > 2, that T~^1Zq C IZ^ and 

7^0 C 7^' C 37^o. 



From Lemma m we have, for 6 eTZ', 



\ A^A^U<r^ 1<r"^■^^<^^ A<rA^U<r^ / 



l<j<k<n l<j<fc<^ l<j<k<n 

where A2 = iA(l + A), ^3 = |A(1 + A)(l + 2A), = ^A(l + A)(l + 6A + 6A 



2^ 



Using Lemma [3l we perform the transformation 6 = T{z). This diagonahzes the 
quadratic form in F{6) and In> becomes: 



1 / 271 



n/2 



F{T{z))dz, 



where 

F{T{z)) =exp(-A252(l + A)-V2 - ^^3^3(1 + A)-V3 - ^^351,2(1 + A)"Vi/^2 
- iA3Si,i,i(l + A)-V? + ^4^4(1 + A)-V4 + ^4^2,2(1 + A)-V2 

+ A4Si,3(l + A)-Vl/i3 - ASi,i,2(l + A)-V?/i2 

+ A4i?i,i,i,i(l + A)-Vt + 0(n-i/2)^^ 

where 



^2 


= - 2, 




i?3 


= n - 4, 




S4 


= n - 8, 




-Bl,2 


= 3(1 -2c) + 


12c 

n 


-Bl,3 


= 4(1 - 2c) + 


32c 

n 


B2,2 


= 3, 






-6c + 12c2 


- 4c' 


n 





-4 + 472 + 0(^-1) 



=^^^ + 0(n'^), 

24c(l - c) ^ 48c2 



5l,l,2 = ^ + ^ = 0(71-'), 

8c^(l-c)(3-c) 8c3(4-c) 
-Di, 1,1,1 = 5 ^ ^ = ^[n 



In order to apply Theorem [2] we choose e = £: + log2/ logn, e' = -^e, e" = e, e'" = 3e 
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e = 0,e = e,A = l, 6{n) = 0(^-^2) ^nd 



(1 + A)2n 2(1 + A) V n 

(1 + A)3n '6(1 + A)2 ^' 
(7- f _ A(l + 2A)(l-v/2) 

^- . -43i?i,i,r^ _ A(l + 2A)(l-v/2) 

(1 + A)3 " ' 3(1 + A)2 ^ 

^ A,B, A(1 + 6A + 6A2) 
^=(rTA)%= 24(1 + A)3 +^^" 



^4^2,2 A(l + 6A + 6A 



(1 + A)4 8(1 + A)3 ' 
(l + A)4ni/2 ^' 

(1 + A)4 ^^"^ ^' 
(1 + A)4 ^^"^ ^' 

^ = [ lem 

2A2 + 2A-1 ~, 

©1 = + 0(n M. 

12A(1 + A) ^ ' 

Theorem now follows from Theorem [21 



4 Concentration of the integral 

In the previous section we proved that the contribution to I{n) from the box TZ' 
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We now consider the contribution to I{n) from the region TZ'^ (defined in (|4])) and show, 
provided A is not too large, that it is neghgible compared to I-jiiin). 

First we import from [5] some useful lemmas. 
Lemma 6. The absolute value of the integrand F{6) of ([3]) is 

\m\= n /(^.+^'^)' 

l<j<k<n 

where 

fiz) = {l + AA2il- cos z)y'^\ 
Moreover, for all real z with \z\ < + X)~^, 

< f{z) < exp{-A2z' + {^A2 + Al)z^). (8) 

Lemma 7. Definet = ^(l + A)^^ and g{x) = —A2x'^ + {jA2+9A2)x'^. Then, uniformly 
for X> and K > 1, 

exp{Kg{x)) dx < ^njiA^K) exp(0(ir-^ + {A^K)-^)). 

-It 

Theorem 8. Suppose that the conditions of Theorem U\ hold, and in addition that 
A = Then 

[ \F{e)\de = 0{n-^)In'{n). 

Proof. The proof follows a similar pattern to that of [HI Theorem 1]. Define t and 
g{z) as in Lemma [71 

Define tt-q, ni, n2, "^,3, functions of 0, to be the number of indices j such that 9j lies 
in [— t, t], (t, vr — t), [tt — t, vr + t], and (— vr + t, — t), respectively. Let TZ" be the set 
of all 6 such that max{non2, C^^), C^^)} > n^"*"^. Any 6 G 71" has the property that 
f{6j + 6k) < /(2t) for at least ra^^^ pairs j, k. Since f{z) < 1 for all z, and the volume 
of TZ" is less than (27r)", we have 

f \F{0)\dO <{2-KYf{2tY'*\ 
Jtz" 

Applying ([8]) and the assumption that A = 0{n'^^^^), we find that 

[ \F{e)\de = 0{e-'^'''^''')In'{n) (9) 

for some Ci > 0. 

For E 1Z'^\1Z" we must have ni, = 0{n^^'^^'^) and either uq = 0(72^/^"^^) or ^2 = 
0(n^/2+e)_ The latter two cases are equivalent, so we will assume that n2 = 0{n^''^~^'^), 
which implies that uq = n — 0{n^^'^^^). 
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Define Sq, Si, S2, functions of 6, as follows. 



m<t}, 

t < \9j\ < 2t}, 
\9j\>2t}. 



Define Sj = 15*41 for each i. Since sq = Uq, we know that Si + S2 = 0{n^^'^~^^). Now we 
bound \F{0) \ in \ TZ" as follows 



fit) < exp 



f{9j + 9k) < I 



A 



if j eSo,ke S2, 



14400(1 + A) 

exp{-A2{9j + 9kf + (^A + Al){9j + 9^)^) if j, k e So, 
1 otherwise. 



Let /2('52) be the contribution to I{n) from those e VJ^yiZ" with the given value of S2, 
and let 0' denote the vector {9j)j^So- The set 5*2 can be chosen in at most n*^ ways. 
Applying the bounds above, and allowing (271)*^"*"*^ for integration over 9j e 5*1 U 5*2, 
we find 

/(S2) <n^'(27r)^^+^^expf 



14400(1 + A) 



I' (so). 



(10) 



where 



H f{9j + 9k)d0' 



j,keSa,j<k 



<[■■■[ exp(-A2 Yl (Oj + Ok)' + (^2^2 + Al) J2 {Oj + Ok)')d0' 
"^"^ j,keSo,j<k j,keSo,j<k 

f- ■ ■ f cxp(-A2(so -2)Y, + 8(so - 1)(^A2 + Al) J2 d0' 

(^j exp(-(so-2)<7(^))rf^) for |(so-l) < 1(^0-2) 
exp(0(l + A-^)n-i) 



< 



< 



< 



< 



^(so-2) 



TT 

A2n 



n/2 



exp 



The third hne of the above follows from the bounds 

p p 
J2 {xj + Xk)'>{p-2)J2^^j and ^ (x,- + x^)^ < 8(p - 1) ^ 

l<j<k<p j=l l<j<k<p j=l 
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valid for aA\ xi,X2, ■ ■ ■ ,xi. The fifth hne follows from Lemma [6l and the last line follows 
from So = n — 0(n^/^+^). Substituting this bound into (1T0|) we find that 

J2I,{s2) = 0{exp{-C2n/logn))I^^in) (11) 

S2>1 

for some C2 > 0. 

With the cases of ([9]) and (iTTl) excluded, we are left with the problem of bounding 
the contribution of G [-2t,2t]" \ 7^^ Let u = 71-^/^+" /{I + A). First note that, by 
arguing as above, we have that 

n 

\F{e)\<l[exp{{n-2)g{e,)) 
i=i 

for all e [-2t,2t]". Also note that, by Lemma H 

2t 

exp((n - 2)g{z)) dz < ^/^^/{A^ exp{0{l + A"^)n"^). (12) 

-2t 

The function g{z) has at most one minimum in [u, 2t], and g(2t) < g{u) for sufficiently 
large n, so 

(^J + J ^ew{{n-2)g{z))dz <Atexp{{n-2)g{u)) 

<exp( , ^ ^ nA. (13) 

- 4(i+A) ; ^ ^ 

Let Ji, J2 be the right sides of (fT2l) and (fT3l) . respectively. Then 

/■ iF(e)ic/0<x:f''V2^""'' 

= jr((i + j2M)"-i) 

= 0(e-^^"^)/7^.(r^) (14) 

for some C3 > 0. 

The lemma now follows from ([9]), ffTTl) and ([H]). □ 



5 Proof of Theorem [T] 

In the case that A = n'~'^^\ Theorem [T] follows from Theorem O and Theorem [HI For 
larger A, the method used in the proof of Theorem [H] is insufficient so we need a new 
approach. 
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Let us assume that we already proved Theorem [T] for A = O(n^). Now we want to 
show that it must be true for larger A as well. First note that for such large A (indeed 
for A/n — )■ oo), Theorem [T] is equivalent to 

Min, i) = V2{X + \r'-'"' (1 + 0(nVA^ + n^)) . (15) 

Let Vn be the polytope of symmetric n x n real non-negative matrices with zero 
diagonal whose rows sum to 1. Then M{n,i) is the number of integer points in iVn- 
That is, M{n,i) is the Ehrhart quasipolynomial of Vn- 

According to Theorem 8.2.6], the vertices of Vn correspond to graphs whose 
components are either isolated edges of weight 1 or odd cycles with edges of weight |. 
That is, the coordinates of the vertices are multiples of ^. By a result of Ehrhart (see 
[21 Ex. 3.25]), there is a polynomial fn{z) with non-negative integer coefficients such 
that 

E^'^(».^)-' = 7r#&T. (16) 

where d is the dimension of Vn- By [Tl], d = n{n—3)/2. 

By applying the binomial expansion to {1 — z'^)^'^^^ in ( IT6|) . we find that M(n, i) is a 
polynomial in i for even i and a possibly-different polynomial in i for odd i. Explicitly, 
there are non- negative integers ho, . . . ,hd (dependent on n and the parity of t) such 
that 

d 



Arguing as in [5], we infer that there is a function a{n,i) such that 

M(n, i) = Q (^Yl (1 + "(^' ^)/^) (17) 

a{nj) > ioi i>d (18) 
a{n, i + 2) < a{n, i) ior i > d. (19) 

Equations ( |T5|) and ( ITTl) both apply for ^ = 6(?2^), so for m G {n^,n^ + 1}, 

MpSn^ = 3^ l + (1 + o(n-'))= 3^(l + 0(n-')). 

M{n,m) 1 + a{n,m)/m^ ^ v v \ J/y 

where the first estimate comes from ( |T7I) and the second comes from ( |T5i) . Comparing 
these two estimates, and noting from (ITSj) and ( |T9|) that < Q;(n, 3m) < a(n, m), we 
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conclude that a{n,m) = 0{n^ ^). By ( IT9|) . this imphes that a{n,£) = 0{n^ ^) for all 
i > n^. Now we can see from (fT7j) that 

M(n, i) = M{n, m) ^}'^'^\ = M{n, m) ( -] (l + 0{n-')) 
M[n, m) \m J 

and apply (fTS!) to M{n,m). This shows that (fTSjl holds for all £ > n^. The proof of 
Theorem [U is now complete. 



6 Naive Thinking 

Define Qx to be the geometric distribution with mean A. That is, for a random variable 
X distributed according to Qx, "we have 

Prob(X=j) = -l-f-A^Y (20) 



1 + A Vl + A, 
for j > 0. 

Define S = S{n,i) to be the probability space of n x n non-negative symmet- 
ric integer matrices with zero diagonal, where each element of the upper triangle is 
independently chosen from Qx- Define events on S: 

Ej : row j has sum i 

n 

E^n ■ Pi Ej 

Eq : the whole matrix has sum n£. 
Note that E^ii C Eq. Also note that each matrix in Eq has the same probability, namely 

1 \©/ A V'/' 

Po ' 



(Proof: Apply fl20|) to each entry in the upper triangle and use the assumed indepen- 
dence of the entries there. The result is independent of the actual matrix entries.) 
Therefore, 

Prob(Eaii) 



M{nJ) 



Pn 



Now make a naive assumption that the events Ej are independent. By symmetry, 
Prob(-Ej) is independent of j, so we get a naive estimate of M{n,i): 

M_e(n,£) = ^^^^. (21) 
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Now consider Prob(£'i). The number of possible first rows is 



(Tliis is the number of ways of writing C. as the sum of n — 1 non-negative integers.) In 
space 5, each such first row has probabihty 



1 V-V A 



Vl+A 

Therefore, 

Substituting this value into ([21]), we get 

Therefore, the formula in Theorem [T] can be written 

M(n,£) = M,,i.e(n,£)v^exp(| + 0(r2-^)). (22) 

Note that v^e^/^ ^ 2.9939. 

After observing a large number of exact values, computed as described in the next 
section, we have noted that fl22|) appears to have an accuracy much wider than proved 
in Theorem [TJ We can even guess extra terms. We express our observations in the 
following conjecture. 

Conjecture 9. For even nl, define A(r;,,£) hy 

M{n, f ) = M^^Un. I) V2 exp + ^ + \ + ^^""'^ 



4 12£(n - 1) n{n - 1) 



Then |A(n,£)| < 1 /or n > 5, £ > 1. 



7 Exact values 

As noted in Section [5l M{n,i) is the number of integer points in iVn, where Vn is 
the polytope defined in that section. Lattice point enumeration techniques such as 
the algorithm in [H] therefore allow the exact computation of M{n, i) for small n. In 
practice this is feasible for n < 9 or with difficulty n < 10, almost irrespective of i. 
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For larger n, Vn has too many vertices for this method to be useful, but we can use 
the technique of [8] and [5]. Define f{z) = 1 + z + z"^ + ■ ■ ■ + z^. Then M{n, i) is the 
coefficient of x{xl ■ ■ -x^i/"^/^ in ni<,<fc<n /(^i^fc?/)- 



If q is any integer greater than max{n£/2, n^(n — l)/2 — ni/2}, then M{n, €) is the 
coefficient of the only term in 



l<j<k<n 



in which each xj appears with a power that is a multiple of £ + 1 and y appears with a 
power that is a multiple of q. Now let p be a prime number for which p — 1 is a multiple 
of both i + 1 and q. Let a and /5 be a primitive {£ + l)-th root and a primitive q-th 
root of unity in Zp, respectively. Then, modulo p, 



q{i + iy 



X 



roH \-rd=n i=0 



^7. • 



9-1 



fc=0 



i=0 



o<«<j<rf 



where the first summation is over all non-negative integers ro, ri, . . . , which sum to n. 
Using sufficiently many primes p, we can extract the exact value of M{n, €) using the 
Chinese Remainder Theorem. As an example of a value computed using this method, 
we have 

M(19, 10) = 613329062511931789477677176839174642138032757885191693120. 

By interpolating in the computed values, we obtain the Ehrhart quasipolynomial for 
small n. Recall that M{n,i) is a polynomial Me{n,i) for even i and a polynomial 
Mo{n, i) for odd i. We have Mo{n, = if n is odd, and the following. 



Me(3,£; 

Me(4, 

M,(5, 

Me(6, 

Mo(6, 
Me(7, 



Mo(4,£) 



256 
19 



120960 



128 
p 9 19 



192 



32 



24 



5376 
857 



I i43_ /7 I ^ «6 I 4567 I 785 «4 , 10919 «3 
4032 24 5760 384 3024 



I 955 p2 I oo_(_ « I 1 
^ 224 ^ ^ 280 ^ 



Me(6,£ 
533 



5 
256 

/)14 



3633315840 



+ 



533 



,13 



86507520 



279413 
2335703040 ' 



fll2 



+ 



9233 
6488064 



on 



3076459 £io _^ 151339 ^9 _^ 4679131 ^8 _^ 9367 ^7 _^ 43502617 



265420800 



2211840 



15482880 



9216 



_^ 478009 £5 _^ 71076539 ^4 _^ 661673 ^3 _^ 1712147 f _^ 



16588800 
9649 



92160 



9123840 



76032 
15 



246400 



2640 



+ 1 



M (R f \ ^ 70241 p20 , 70241 pi9 , 18703309 pis 

ivie{o,f,) 5088422500761600 "^ 72691750010880 "^ 585359881666560 

, 12330581 pi7 , 428460619 pid , 33009749 pi5 
"^ 18582853386240 "^ 44144787456000 "^ 310542336000 

, 90842880341 nU , 5580172163 /,13 , 1110632463421 fi2 
"^ 100429391462400 910924185600 33108590592000 

, 4381892419 nil , 304644862903 «io , 22001378209 «9 , 262880239845943 /i8 
"^ 29196288000 551809843200 13138329600 62768369664000 

, 12867890603299 «7 , 3890196991231 p6 , 9530810537 n5 , 76295531167 /)4 
"^ 1494484992000 ^ 269007298560 485222400 "^ 3592512000 

, 1100694281 -73 , 50787821048 f2 , 135038369 f , -. 
"^ 61776000 "^ 4583103525 "^ 29099070 "^ 

Mo{8,i)^Mei8,e) 

35 /i5 1225 M 13685 /,3 17885 /i2 233261 « 78057 



1048576 2097152 3145728 1048576 6291456 2097152 

71^ /q /^ ^ 863924282670630091 f27 

' ^ 7732694804887618394297204736000000 

_^ 863924282670630091 ^^26 , 10311705659720524879 ^25 



+ 



71599025971181651799048192000000 ' 16522852147195765799780352000000 

44159888290330963 ^24 , 44603828594214317123 ^23 

2145824954181268285685760000 91793623039976476665446400000 



, 4134171051301720697 f22 . 2139768518991928638127 f2i 

"^472621628924364010291200000 17143275449165567282380800000 
, 2365877475528196499 p20 , 167364777037473990001 ^19 



1632692899920530217369600 ' 12005094852356839833600000 

43210221809651966023 «18 , 7598908879241416557943 «17 
383621452048996761600000 ^ 9846283935924250214400000 

78473046995519797477 /,16 , 2690417378247820105229333 pi5 
17375795181042794496000 ^ 118589802110617072435200000 
12598164604216578106061 /)14 , 39802237244716247322233 /<13 



+ 
+ 

"^ 128343941678157004800000" ' 108598719881517465600000" 
, 183315648883655207683 ni2 , 11492891877126624163867 nil 
"^ 155141028402167808000 3496549693154918400000 
, 20646561932994651460327 /-lO , 12699041960623534314853039 «9 
"^ 2622412269866188800000 "^ 784756871757456998400000 
, 3536936635157608410019 fS , 602776622158017864239297 n? 
"^ 124564582818643968000 14273025114636288000000 
, 8959111748174759872739 /<6 , 62149609860286754066479 f5 
"^ 169916965650432000000 "^ 1139859644571648000000 
, 416558573311485749 /,4 , 7739053944610908107 «3 
"^ 9089789829120000 254233401117696000 
, 1309315468639693 ^2 , 94565099767 f , -, 
"^ 85753329742080 "^ 17847429600 "^ 

The same method would yield M{10,£) with a considerable amount of computation. 
For completeness, we also give the Ehrhart series Ln{x) = ^£>o M{n, t)x^ for n < 9. 

(1 - x')L^{x) = 1 
(1 - xfLi{x) = 1 

(1 - a:2)^L5(a;) = {x^ + 1) + 16 {x^ + x^) + 41 a;^ 
(1 - x)^°(l + x)L6(x) = (x^ + 1) + 6 (x^ + x) + 30 (x^ + x^) + 40x2 
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(1 - x^Y'^Lrix) = + 1) + 807 (x^^ + x^) + 81483 (x^^ + x^) 
+ 1906342 (x^^ + x^) + 15277449 {x^^ + x^) 
+ 50349627 (a;^^ + + 74301542 

(1 - xy\l + xfLsix) = {x^° + 1) + 90 {x^^ + x)+ 4726 {x^^ + x^) 

+ 107050 (x^^ + x^) + 1261121 (x^^ + x^) 
+ 8761248 (x^^ + x^) + 39187016 (x^^ + x^) 
+ 119662536 (x^^ + x^) + 259344246 (x^^ + x^) 
+ 408811676 (x^^ + x^) + 475095180 x^° 

(1 - x^Lgix) = (x^^ + 1) + 52524 (x^^ + x') + 169345602 (x^^ + x^) 
+ 78276428212 (x^^ + ^ 10217460516057 (x^° + x^) 
+ 527531262668208 (x^^ + x^°) 
+ 13016462628712186 (x^^ + x^^) 
+ 172410423955058664 (x^^ + x^^) 
+ 1322251960254170931 (x^^ + x^^) 
+ 6176715510750440488 (x^° + x^^) 
+ 18182086106689738044 (x^^ + x^O) 
+ 34470475812807166836 (x^^ + x^^) 
+ 42606701216240491693 x^^ 

Machine-readable versions of these exact formulas, along with many other exact 
values of M{n,£), can be found at |10] . 



8 Concluding remarks 

In this paper we have begun the asymptotic enumeration of symmetric non-negative 
integer matrices with given row sums, by considering the special case of uniform row 
sums and zero diagonal. Further cases, which can be approached by the same method, 
are to allow the row sums to vary, and to allow diagonals other than zero. The structure 
of random matrices in the class can also be investigated by specifying some forced 
matrix entries. We hope to return to these problems in the future. 
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